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Abstract 

This analysis aims at exploring what can be said about the growth rate of magnetized in- 
homogeneities under two concurrent hypotheses: a phase of quasi-de Sitter dynamics driven 
by a single inflaton field and the simultaneous presence of a spectator field coupled to gravity 
and to the gauge sector. Instead of invoking ad hoc correlations between the various compo- 
nents, the system of scalar inhomogeneities is diagonalized in terms of two gauge-invariant 
quasi-normal modes whose weighted sum gives the curvature perturbations on comoving or- 
thogonal hypersurfaces. The predominance of the conventional adiabatic scalar mode implies 
that the growth rate of magnetized inhomogeneities must not exceed 2.2 in Hubble units 
if the conventional inflationary phase is to last about 70 efolds and for a range of slow roll 
parameters between 0.1 and 0.001. Longer and shorter durations of the quasi-de Sitter stage 
lead, respectively, either to tighter or to looser bounds which are anyway more constrain- 
ing than the standard backreaction demands imposed on the gauge sector. Since a critical 
growth rate of order 2 leads to a quasi-flat magnetic energy spectrum, the upper bounds on 
the growth rate imply a lower bound on the magnetic spectral index. The advantages of the 
uniform curvature gauge are emphasized and specifically exploited throughout the treatment 
of the multicomponent system characterizing this class of problems. 
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1 Formulation of the problem 



According to a recurrent theme of speculations, large-scale magnetic fields could be generated 
in the early Universe [H El [3l IH El E] . The curvature perturbations evolving for typical length- 
scales larger than the Hubble radius can thus be magnetized with a mechanism bearing some 
resemblance to a pristine non-adiabatic pressure fluctuation. This observation has been used 
some time ago to argue that the evolution of curvature perturbations constrains the magnetic 
power spectra [7]. In the present paper the same logic explored in [7j, i.e. the predominance 
of the adiabatic mode over the gauge contributions, will be used to analyze consistently the 
fluctuations of inflationary magnetogenesis and derive different constraints on the growth 
rate of the corresponding inhomogeneities. 

The fate of magnetized scalar modes during diverse dynamical regimes can be followed 
through a gauge- invariant variable, conventionally denoted by (, describing either the curva- 
ture perturbations on the hypersurface where the energy density is uniform or, complemen- 
tarily, the density contrast on uniform curvature hyper surf aces. The latter interpretation 
becomes physically appealing and mathematically simpler in the so-called uniform curva- 
ture gauge which has been discussed in different contexts [H [9l [TOl [H] . Since ( is ultimately 
gauge- invariant its evolution can be studied in any gauge and the result of [7] , derived orig- 
inally in the uniform curvature gauge [12 , can be confirmed in different coordinate systems 
and in different dynamical situations [13] □; neglecting electric fields and Ohmic currents the 
evolution equation of ( is: 

C = -^^^Pnad + ^^'f~% B - % (1.1) 

Pt{l + Wt) 3pt{l + wt) 3 

where Spna.d{x,T) accounts for the non-adiabatic pressure inhomogeneities; 6pb{x,t) is the 
fluctuation of the magnetic energy density and 9t = V ■ Vt is the divergence of the total 
velocity field. Barring for a possible contribution of the total velocity fielclf] and in the 
absence of entropic modes (i.e. 5pnad = 0) the solution of Eq. (II. ip is in fact a functional 
of the total barotropic index Wt = Pt/Pt and of the total sound speed c^^ = p'Jp\. Denoting 
with C*(2^) the conventional adiabatic mode, the full solution of Eq. ( II. ip becomes: 

C(x,a,a*) = C*(x) + / . —bpB{x,h)d\iih, (1.2) 

J a, 3pt(o)[l + Wt(o)J 



where the integration variable is provided directly by the scale facto 



^As usual the prime denotes the derivation with respect to the conformal time coordinate t and "H. = a! j a 
where a is the scale factor of a conformally flat metric of Friedmann-Robertson- Walker type. 

•^This term is subleading for wavelengths larger than the Hubble radius at the corresponding epoch. 

"^Equation corresponds exactly to Eq. (2.15) of Ref. [7j. The same equation has been used [H] 

to deduce the initial conditions of the Cosmic Microwave Background (CMB) anisotropics in the presence 
of postinflationary magnetic fields characterizing the so called magnetized adiabatic mode. Equation (jl.ll) 
has been later generalized to the case when Ohmic currents are present and also in the framework of the 
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In the uniform curvature gauge [SI El [IDl [H] , Eq- fll-ip stems directly from the covariant 
conservation of the total energy-momentum tensor on uniform curvature hyper surf aces: 



5p[ + m{5p, + 5pt) + (pt + Pt)^t = 0, 



c 



{5pt + 5pb) 
3(1 + wt)pt 



(1.3) 



where, by definitioijj, 5pt = cl^5pt + (5pnad- The covariant conservation of the total energy- 
momentum tensor also implies the adiabatic suppression of 5p-b{x^t) redshifting as a~^. A 
more general derivation of Eq. fll.ip including Ohmic currents and energy flow is swiftly 
outlined in Eq. (1A.28P of appendix \K\ 

The strategy leading to Eqs. (Il.ip - (ll.3p could be rigidly translated, at first sight, directly 
during the inflationary stage of expansion. It might then seem plausible to keep the whole 
logic untouched but to concoct specific modifications of the evolution of dpB modeling, via an 
appropriate rate of increase, the growth of 5p-Q during inflation when the relevant wavelengths 
of the corresponding fluctuations are larger than the Hubble radius. A candidate equation 
describing the amplification of the magnetic inhomogeneities is, for instance. 



where 2J^ denotes the rate of increase of the magnetic energy density which is twice the 
growth rate of the magnetic field itself. Barring for the presence of Ohmic currents and 
electric fields, Eq. (11. 4p partially accounts for the effect of superadiabatic amplification of 
the magnetic fields but disrupts the covariant conservation of the total system. This means 
that the evolution equation for C, is no longer valid. A compensating term can be added 
at the right hand side of Eq. (II. 3p but this has different drawbacks since the evolution 
equations derived from the covariant conservation of the total energy-momentum tensor will 
be no longer compatible with the remaining perturbed Einstein equations. 

If the dynamics of the inflationary magnetogenesis is not taken into account specifically, 
the evolution of the whole system turns out to be inconsistent because of the lack of covariant 
conservation of the total energy-momentum tensor. The first mandatory step for any analysis 
involving the fluctuations of inflationary magnetogenesis is to posit a perfectible framework 
where magnetic fields are amplified, the Bianchi identities are satisfied and the inflation- 
ary dynamics is satisfactorily implemented. We suggest that the dynamics of magnetized 
inhomogeneities can be consistently scrutinized in the following system: 



gradient expansion (see, respectively, the first and second paper of [13]). Exactly the same equation 

has been applied in Ref. |14) with virtually the same purpose of deriving a bound connecting the amplitude 

of the adiabatic mode and the strength of the magnetic field. 

^The total pressure can fluctuate either because of a change in the energy density (when the specific 
entropy is unperturbed) or because of a change in the specific entropy of the system (when the energy 
density is unperturbed). 



5p'b + AUSpB = 2T5pB, 




= snG t;{^) + T;(a) + r;(p, p) + z;(y) , 
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denotes the Einstein tensor while T^(v5) is the energy-momentum tensor of the inflaton 
(p; Tl^{cr) is the energy- momentum tensor of a spectator field a and T^^ip^p) is the energy- 
momentum tensors of the total fluid sources while Z'j^^iY) is the energy- momentum tensor of 
the gauge fields. The explicit coupling to the spectator or to the inflaton fields leads to the 
covariant non-conservation of Z]", 

= ^ Yo^p Y^^ + f y,,, (1.6) 

where Yap is the gauge field strength, ja is the four-current and A (a;) is a function parametriz- 
ing the coupling between the gauge fields and the spectator field a. For sake of generality 
we shall also consider the possibility the coupling will depend both on a and so that 
A = A((7, The covariant non-conservation of Z^^ is compensated by the covariant non- 
conservation of the other energy-momentum tensors: 

V;,7;" = -j<'y„„. (1.9) 

Equation (11. 5p captures a class of magnetogenesis scenarios studied along different perspec- 
tives through the years and some of the possibilities will now be recalled. In general terms 
A = X[ip{x), cr(3;), ...] niay be a functional of various scalar degrees of freedom such as the 
inflaton ip [15] , the dilaton [TH] , a dynamic gauge couphng [T71 HB] (see also [121 ED] ) • The 
field A can be a functional of a spectator field cr, [211 122] (see also [23l [21]) evolving during 
the inflationary phase; in this case there is no connection between the evolution of A and the 
gauge coupling. Some of these possibilities can be realized in the case of bouncing models 
[16] , some other are compatible with the standard inflationary paradigm [151 [13 [HI 121] • It 
is finally worth recalling a recent observation: the initial conditions of inflationary magneto- 
genesis may be conducting [27] since the Ohmic currents present during the preinflationary 
dynamics are not damped by expansion due to the Weyl invariance of the electromagnetic 
sources. 

A perturbative treatment of the fluctuations of inflationary magnetogenesis in a consistent 
dynamical framework encompassing the inhomogeneities of the inflaton, of the spectator 
field, of the growth factor and, last but not least, of the relevant plasma variables will now 
be presented. This analysis is lacking and it is mandatory if the principle of predominance of 
the adiabatic mode, spelled out in of [7], is to be enforced during the inflationary phase. The 
tools developed in this paper will allow for an accurate constraint involving simultaneously 
the slow roll parameters, the total number of inflationary efolds and the total rate of increase 
which can be defined, for the present purposes, directly from Eqs. (II. 6p as J-" = drVX/V^- 
Inspired by the analysis of [7] the pivotal variables for the evolution of the gauge sector will 
not be the gauge fields but rather the components of the energy-momentum tensor. This 
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strategy together with the gauge choice mentioned above will allow for a swifter calculation 
of the primary and secondary curvature perturbations induced by the infiaton field and by 
the spectator field. 

The layout of this paper is the following. In Sec. |2]the general equations of the system will 
be discussed and the main notations specified. In Sec. Othe description of stochastic averages 
will be introduced with the aim of reducing the evolution of the system to the evolution of the 
components of the energy-momentum tensor of the gauge field inhomogeneities. In Sec. H] 
the quasi-normal modes of inflationary magnetogenesis will be discussed in general terms. In 
Sec. [5]the magnetized power spectra of the scalar modes will be computed while in Sec. [6]the 
bounds on the growth rate of the magnetic energy density will be derived. To avoid lengthy 
digressions various technical details have been collected in the appendices: in appendix |X] 
the evolution equations of the system have been explicitly derived in the uniform curvature 
gauge systematically used in the analysis; in appendix |B] the second-order correlations of the 
electric and magnetic fields have been specifically computed and analyzed. 



2 Basic equations and definitions 



The system of equations swiftly outlined in Eqs. (11. 6p and (ll.7p - (ll.9p can be illustrated by 
specifying the actions of the different contributions: 



Scy 



tot '-'gravity ~l~ -\- Sfj -\- S^^a ~l~ 'S'fluid) 



where the first three terms of Eq. (12.11) are given by: 



So, 



gravity ~l~ S^ 



2P'~ ' 2 



^jd'x^ A(v9, a) V y'"" - I d'x^j, + + S^^y 



(2.1) 

(2.2) 

(2.3) 
(2.4) 



In Eqs. (I2.2p -( |23|) . V{ip) and W{a) denote, respectively, the potentials of the infiaton field 
and of the spectator field. In Eq. (12. 4p = jj^^ — is the total current; S(±) are the 
actions of the charged species while the last term of Eq. (12.11) (parametrized via a barotropic 
fluid) can be important either at the onset of inflation (for conducting initial conditions |27j ) 
or during the postinfiationary phase. The notations for the Planck length and for the Planck 
mass in units ^ = c = kb = 1 are as follows 



= 8 TT G 



»7r 
Ml 



M 



■2 ' 



(2.5) 



where Mp = G'^/^ ^ x 10^^ GeV. On top of Eq. ([L5]), The equations of motion of the 
various fields appearing in Eqs. (12. ip and (12. 4p are given by Eq. ( 11. 5p supplemented by the 
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following three equations: 



dV 



1 dX. 



da IQtt da 



V^Tp^ = 0; 



the explicit forms of the energy-momentum tensors T^{(p), T^{a) and T£{p, p) are: 



T^{a) = d^ad^a- 



^g^'d^ad^a - W{a) 



5^ 



47r 



(2.6) 

(2.7) 
(2.8) 



(2.9) 

(2.10) 
(2.11) 
(2.12) 



where g^^ UaUp = 1. By using Eqs. fl2.9p - fl2.10l) and fl2.lip - p.12p the evolution equations 
for the energy-momentum tensors mentioned in Eqs. fll. 70 - 01. 9p . Eqs. (12. 6 P and (12. 7p can 
be reproduced bearing in mind the following pair of equations for the gauge fields: 



0/3 



0; 



(2.13) 



the dual field strength is defined as Y""^ = E°''^'^'^Y^^/2 in terms of the Levi-Civita tensor 
density E°''^'^^ = e"^^^ / \/—g- Note, finally, as already mentioned in Sec. [Uthat A = A(<y9,cr) 
and, consequently, 9^A = {d^Xd^ip + da\dfj,a). 



2.1 Background evolutions and some approximations 

In a conformally flat background of the type g^j^ = a'^{T)r)ai3 (where a(r) is the scale factor 
and ?7q,/3 is the Minkowski metric), Eqs. (11. 5p and fl2.6p - fl2.7p lead to the following set of 
equations valid during the inflationary phase 



2Mp {n^ - W) = <^'' + a'\ 
dV 

if" + 2n^' + — a' = 0, 

dip 

a" + 2na' + ^a' = 0. 
da 



(2.14) 
(2.15) 
(2.16) 

(2.17) 



As mentioned prior to Eq. (11. ip . in Eqs. fl2.14p - fl2.17p the prime denotes a derivation with 
respect to the conformal time coordinate r; furthermore T-t = (In a)' = aH where H = a/a 
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is the conventional Hubble rate and the overdot denotes a derivation with respect to the 
cosmic time coordinate t. 

The slow roll approximation completely defines the evolution during the infiationary 
phase where the parameters e, t] and fj are all much smaller than 1 and eventually get to 1 
when inflation ends. The definitions of the slow roll parameters within the notations of this 
paper are as follows: 



note that and V<p<p are shorthand notations for the first and second derivatives of the 
potential V{(f) with respect to (p. The slow roll parameters 77, r] and e are not independent 
and their mutual relation, i.e. t] = e — r], follows from the slow roll equations written in the 
cosmic time coordinate: 

f)Y , „ . 

3H^ + — = 0, 3MpH^ = V, 2MpH = -0^ (2.19) 

where, by definition of spectator field, we have that po- ^ and 0^ ^ having intro- 
duced the energy densities of the infiaton p^^ and of the spectator field Po-- In the slow roll 
approximation and for constant e we have that 

'H = aH = - . \ . (2.20) 
(1 - e)r 

There are some classes of exact solutions which shall be used in order to test the specific 
approximations discussed in the second part of this analysis. If both (p and a have exponential 
potentials a solution of the system f l2.14p -f l2.17p subjected to the constraint that p^- ^ Pip 
and 0^ ^ (j^ can be written, in cosmic time, as: 



a{t) = {Hi ty, ip{t) = V2a Mp In [Hi t), (2.21) 



Viip) = M;Ht{3a'-a)exp 
a(t) = 2Mln(Mt), iy((T) = 2(3a - 1)M^ exp 



(2.22) 
(2.23) 



with M ^ Mp and a ^ 1 so that e = r] = 1/a 1. In conformal time the corresponding 
scale factor becomes: 

a(r) /3 = ^, (2.24) 

V Ti / a — 1 

with /3 — )■ 1 in the limit a ^ 1 and e ^ 1. In specific models of inflationary evolution, 

the values of the slow roll parameters, for a given number of efolds, can be related to the 

properties of the potential. To keep the discussion sufficiently general we shall treat the slow 

roll parameters and the number of efolds as independent variables; conversely, as already 

mentioned prior to Eq. fl2.20p . the slow roll parameters will be taken to be constant implying 

that the inflationary potentials considered here have a monomial form. 
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3 Quantum and stochastic descriptions 



The amplification of tlie gauge fields can be described quantum mechanically in terms of the 
appropriate canonical field operators and of their related mode functions. This description 
is equivalent to the evolution of the power spectra of the different correlations. The two 
approaches are related and this observation turns out to be very practical for the present 
considerations. 



3.1 Evolution of the canonical gauge field fiuctuations 

In the conformally flat background discussed in the previous section, Eq. fl2.13p becomes 
explicit in terms of the canonical electric and magnetic fields diagonalizing the action and 
the canonical Hamiltonian 



-Lv . (v^^) = 47rgK - v^V • = 0, (3.1) 

-^V X (VXB) =47rq{n+v+ - n^v^) + l=-^{^/XE), (3.2) 

— # — * 

where E(x, r) and B(x, r) are 

E = a^VXe, B = a'^y/\b. (3.4) 

The fields e and b are introduced from the corresponding field strengths, i.e. Yio = —a?ei 
and Yij = —a^eijkb^. The gauge action is canonical in terms of E and B and not in terms 
of e and b. Furthermore the system of Eqs. ( I3.ip - fl3.3p . in the absence of electromagnetic 
sources, is invariant under the generalized duality transformation E — )• —B, B ^ E and 
y/X^l/y/X [251 [26] (see also the second paper quoted in Ref. [27]). 



3.2 Evolution of the power spectra 

Let us start by recalling the notion of stochastically distributed Fourier modes in the case 
of the electric and magnetic fields, i.e. 

{B,iq,T)Bjip,T)) = '^P^(q,r)P,,iq)S^'\'l + P). (3-5) 

2vr2 

{Ei{q,T)E,{p,T)) = —P^{q,r)P,,{q)S'^'\q + p), (3.6) 
where Pij{q) = {Sij — qiqj) (with qi = qi/\q\)', the conventions for the Fourier transform are: 
B^{x, t) = J £k B,{k. r) e-^^-^ E,{x, r) = J d'k E,{k, r) ^^'^ (3.7) 
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The evolution equations fl3.ip . fl3.2p and fl3.3p are equivalent to the following set of equations 
obeyed by the power spectra of Eqs. fl3.5p - fl3.6p 

= 2^Pb - g^EB, (3.8) 

dP 

= -2( J- + 47rae)PE + qPeb, (3.9) 

= 2g(PB - Pe) - 47rorePEB, (3.10) 

where cxc denotes the Ohmic conductivity, J-" = / \/\ is the growth rate and Pee is the 
cross- correlation spectrum defined implicitly by the following equation 

{E -V X B) + {B -V X E) = 2 J dqP^Biq^r). (3.11) 

The cross-correlation spectrum provides the physical difference between a stochastic collec- 
tion of gauge fields (described by Eqs. ( 13. Sp and (13. 6p ) and their quantum analog which 
will be discussed in a moment (see Eqs. ( I3.17p -( 13.17p ). Conducting initial conditions [27] 
correspond, in Eqs. ( I3.6p -( l3l8|l . to the limit Per where the magnetic fields are amplified 
and the electric fields suppressed either at the same rate or even exponentially depending on 
the value of the protoinfiationary conductivity. In quantum mechanical terms the canonical 
normal modes are field operators defined a^ 



(3.12) 
(3.13) 



(27r)^/^ a ^ 

where the evolution of the mode functions is given by: 

fk = J^fk - 9k, 9'k = -^9k - 47r gu + /^^ /fc, (3.14) 

and the possibility of conducting initial conditions has been included for comparison. In the 
absence of sources, as already mentioned after Eqs. (I3.ip - (l3.3p . Eqs. (I3.14p are invariant 
under generalized duahty transformations stipulating that fk — )■ 9k/ Qk — ?• —kfk and 

In terms of the mode functions, the Fourier components of Bi{x,T) and Ei{x,T) are 
respectively 

Bi{q, r) = -iemni ^ qm[a,^,a fq{r) + al^.^^ /*(r)], (3.15) 

a 

m,T) = 9,{r)+al.^ g;{r)]. (3.16) 

/3 



®Note that e\°'\k) (with a = 1, 2) are two mutually orthogonal unit vectors which are also orthogonal to 
k; furthermore J2a = PijW- 
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It can be immediately checked that Eqs. f l3.15p - fl3.16p obey the stochastic averages defined 
earher in Eqs. fl3.5p -f l3l6|) : for instance, in the case of the magnetic field operator, 



(0|5,(g,r)4(p,r)|0) 



g3 



2^2 1^9 



\Ur)\\ (3.17) 



in full analogy with Eq. fl3.5p . It can be easily argued that Eqs. f l3.6p and f l3.8p -f 1XTU|) are 
similarly satisfied with 



^E(g,r) 



27r' 



■\g,{T)\\ PEB{q,r) 



271 



,[f:ir)g,{T) + f,{T)g;{T)l 



(3.18) 



where Pe(9, t) denotes the power spectrum of the electric fields and Peb{(1, t) is the spectrum 
of the cross-correlation between electric and magnetic fields. To have compatibility between 
the evolution equations of the power spectra (i.e. Eqs. f l3.8p -f lXTU]) ) and the evolution 
equations of the mode functions (i.e. Eq. fl3.14p ) the cross-correlation spectrum is essential. 
Using the power spectra defined earlier the average magnetic and electric energy densities 
are 

Backreaction problems are avoided if and are smaller than the background energy den- 
sity 3H^Mp. Moreover the contribution of the electric and magnetic fields to the evolution 
equations of and a must be subleading. These requirements are, however, less severe than 
the ones stemming from the predominance of the adiabatic mode discussed in Sec. [61 



3.3 Fluctuations of the energy- momentum tensor 

The normalized fiuctuation of the energy density are given by 



6pBix,T) 
where 



(27r)3/2 



SpB{q,T) e 



iq-x 



6pe{x,t) 



d^q 



(27r)3/2 



5pE(g,r)e-«^ (3.20) 



(27r)3/2 
1 

(27r)3/2 



d^k 
d^k 



B,ik,r)B,{q-k,T) 
Ei{k,T)E,{q-k,T) 



47r2 
47r2 



PB{k,T)6^^\q) 
PE{k,T)S^'\q) 



(3.21) 



Similarly, the normalized fiuctuations of the electric and magnetic pressures are 6pb{x, t) = 
6pb{x,t)/3 and Spe{x,t) = 6pe{x,t)/3. The electric and magnetic anisotropic stresses are 



nSf(f,r) 



(27r)3/2 



1 



(27r)3/2 



d'qn^\q,T)e- 



■iq-x 



(3.22) 
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where 



nSf(g,r) 



1 



J (27r)3/2 

^^'""^ 7 (277)3/2 



r)5, (g'- k, t) - -fBm{k, r)B^{q- k, r) 
Ei{k, T)Ej{q- k, r) - -fE^ik, T)E^{q- k, r) 



,(3.23) 
.(3.24) 



It is practical to introduce the scalar projections of the electric and magnetic anisotropic 

stresses 

V2nB(x,r) = 9,9,n'j^(x,r), V2nE(x,r) = 9,9,ng(x,r), (3.25) 

entering the evolution equations of the scalar modes of the geometry. Note that the stochastic 
averages of the fluctuations variables defined in Eqs. fl3.20p - fl3.2ip and fl3.22p - fl3.24p are all 
vanishing, i.e. using Eqs. fl3.5p - fl3.6p . {6pb{x,t)) = and {6pE{x,r)) = 0. The second order 
correlations of the energy density fluctuations and of the anisotropic stresses are defined as 



{Spx{q,T) 6px{p,r)) = -^Qx{q,r) 6^'\q + p) 



27r2 



{Ux{q,T)Ux{p,r)) = Q^^{q , r) S^'^Xq + p) , 

q^ 

where X = B, E leading, overall, to four independent spectra 



(3.26) 
(3.27) 



QB{q,T)-- 
QEiq,r)-- 
QBn(g,r) 



128 7r3 

,3 



q- 



128 7r3 

,3 



d^k 



d^k 



PB{k,r)PB{\q-k\,T) 



fc3 



\q — A;|3 



PEik,T)PE{\q-k\,T) 



k^ 



\q — fc|3 



Ap(A;,g), 



^p{k,q), 



288 7r3 a8(r) 

,3 



288 7r3 a« r 



d-^k 



d^k 



PB{k,T)Pj,i\q-k\,T) 



fc3 



\q — k\ 



PE{k,r)PE{\q-k\,T) 



fc3 



\q — k\ 



^u{k,q), 



An{k,q). 



(3.28) 
(3.29) 
(3.30) 
(3.31) 



The functions Ap{k, q) and An(A;, q) are defined as 
Ap{k,q) 

Au{k,q) 



^ I [k-{q-k)] 
k'^\q — fc|2 
1 I [k-{q-k)f ^ 6 

fc2|g*— fc|2 g2 



(3.32) 



k ■ {q — k) 



[k ■ (g - 
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fc2|g-fc|2-2[A?-(g-A?)|]2 + 



A;2|g'— A;|2 



(3.33) 



The functions Ap{k,q) and An(fc,5') coincide for magnetic and electric degrees of freedom 
since both E and B are solenoidal fields: B is solenoidal because of the absence of magnetic 
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monopoles while E is solenoidal because the pprotoinflationary plasma is globally neutral 
and any electric charge asymmetry is absent. The explicit expressions of the power spectra 
of Eqs. fl3.28l) - fl3.31l) are presented and discussed in appendix [B] in the case of a monotonic 
growth rate. 



3.4 Eenergy- momentum tensor evolution 

Instead of computing 6pB, 5pE and the other relevant components from the field variables, 
it is practical to deduce the evolution equations obeyed by these quantities. The various 
components of the energy-momentum tensor can be written, in real space, as 

^0 = PB+PE + SPB + 5PE, (3.34) 

2i = -(j^^+PB}Si-{6p^ + 6pB)6i + uf^' + uf^'. (3.36) 

Using Eqs. fl3.34p -f l3.36p and the explicit expression of the covariant derivative, Eq. (11. 6p 
demands, in components, 

driSpE + Spb) + 4:H{6pE + 6pb) = 27{8pB - 8p^) -P- (3.37) 

d^P + mP = _ v2[5pB + - (Hb + He)], (3.38) 

where P = V ■ S denotes the three-divergence of the Poynting vector S = {E x i?)/(47ra^). 
In Eq. (13.381) the terms di[{d'^ X) / X\{S pB — 5pe) and di\d\5pB — 6pe)/X have been neglected 
since they couple spatial gradients of the growth rate and magnetic inhomogeneities. These 
terms are of higher order in the present description. Furthermore, using standard vector 
identified Eq. (13.381) can be recast in the following form: 

a.p^4np = y-^': a) ^ V- [(V X B) X g| + y. |(v x e) x e]^ 

The evolution of the difference between 6pB and 6pE can be obtained directly from Eqs. 
fl331)-fl33|): 



dr{6pB - 6pe) + m{6pB - 6pE) = 2F{6pE + 6pb) - B'^^E + E-VxB ^ E^ 

The system of Eqs. f l3.37p -f !3.40p can be studied in various approximations (subleading 
spatial gradients, large conductivity limit and so on and so forth). In the most naive case 
P{x, t) simply scales as a~^. This can be easily understood since, up to spatial gradients, 
the evolution of P does not depend on the growth rate. Conversely, the time derivative of 
P is proportional to the Laplacians of the pressures and of the anisotropic stresses. 



^Givcn a solenoidal vector field d, (such as B or E) the product di Cjd^C^ can be expressed as V • [(V x 
C) X C] + V2CV2. 
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4 Quasi-normal modes 



The evolution of the scalar modes of the geometry, of the inflaton and of the spectator field 
are all coupled to the scalar inhomogeneities of the gauge sector. This system will now 
be reduced to the evolution of its quasi- normal modes whose equations are coupled but, 
most importantly, decoupled from all the other perturbation variables. The considerations 
of the present section and of the appendix |X] can be easily generahzed to various situations 
involving, for instance, more than one spectator field. 

4.1 Uniform curvature hypersurfaces 

The scalar fluctuations of the four- dimensional metric are parametrized by four different 
functions whose number can be eventually reduced by specifying (either completely or par- 
tially) the coordinate system: 

^sfl-oo = 2a^0, 6sgij = 2a^{tp6ij - didja), S^goi = -a^dif3, (4.1) 

where 6s denotes the scalar mode of the corresponding tensor component; the full metric 
(i.e. background plus inhomogeneities) is given, in these notations, by gai3{x,T) = gap{T) + 
6sgai3{x,T) where, as already mentioned prior to Eqs. fl2.14l) - fl2.17l) ^q,^(t) = a'^{r)rjaj3- For 
infinitesimal coordinate shifts r — )■ r = r -|- eo and — >■ = -|- d^e the functions 0(x, r), 
I3{x,t), 'ijj{x,T) and a(x, r) introduced in Eq. (14. ip transform a^ : 

0^0 = 0_7^eo-e[„ ^^^ = ^ + neo, (4.2) 
/3-)-^ = /3 + eo-e', a = a - e. (4.3) 

In the uniform curvature gauge two out of the four functions of Eq. fl4.ip are set to zero 

mm- 

a = 0, ?/' = 0, = 0(f,r), i3 = i3{x,t). (4.4) 

Starting from a gauge where a and ip do not vanish, the perturbed line element can always 
be brought in the form (14.40 by demanding a = and -0 = in Eqs. (14. 2 p and (14. 3p . If 
a and ip ^ 0, the uniform curvature gauge condition can be recovered by fixing the 
gauge parameters as e = a and eo = —ip/T-L. This choice guarantees that, in the transformed 
coordinate system, ip = a = 0. 

A convenient gauge choice is essential for a sound treatment of problems involving the 
presence of anisotropic stresses. The conformally Newtonian gauge is known to be unsuitable 
for the analysis of perturbative systems where the anisotropic stresses play an important 
role. Similar caveats arise in the discussion of the Einstein-Boltzmann hierarchy whenever 

^The slow roll parameter e must not be confused with the parameter of the gauge transformation. These 
two variables never appear together either in the preceding or in the following discussion so that no confusion 
is possible. 
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the entropic initial conditions are dominated by the anisotropic stresses as it happens in 
the neutrino sector. Both points have been addressed long ago when discussing the initial 
conditions for the magnetized CMB anisotropies [28] (see also [12] and references therein). 
Other gauges are also suitable for the treatment of magnetized inhomogeneities but we shall 
not discuss them here. As already mentioned in Sec. [H to avoid lengthy digressions the full 
set of evolution equations has been presented and discussed in the appendix El 



4.2 The decoupled system 



Consider, to begin with, the evolution equations for the fluctuations of the infiaton (i.e. x<fi) 
and of the spectator field (i.e. Xa) which are reported in Eqs. flA.20p and flA.2ip . Recalling 
that Act = V^Xcr ^iid A,^ = V^Xv?) from the momentum constraint of Eq. (1A.9P (neglecting 
a generic fluid contribution which is anyway irrelevant during the inflationary phase) the 
following relation holds: 



A. 



-P. 



(4.5) 



During inflation the three-divergence of the Poynting vector P decreases always as so the 
predominant contribution to the curvature perturbations on uniform curvature hypersurfaces 
is given by the first two terms of Eq. (14. 5p . There is, however, an important proviso: the 
time derivative of A^ (i.e. A'^) appearing in Eqs. ( ]A.20p and ( ]A.2ip leads to a term going as 
P' = drP containing the Laplacians of the magnetic and electric energy density fluctuations 
(see Eq. fl3.38p ). It is advisable, as usual, to assess the relative weight of different terms not 
at the beginning, but rather at the end of the derivation. 

In the gauge (14. 4p . the curvature perturbations on comoving orthogonal hypersurfaces, 
customarily denoted by IZ, coincides with up to a background dependent coefficients. 



7^ 



Defining A-ji = V^TZ and recalling Eq. (14. 6p . we have, from Eq. (14. 5p . 



A 



Ha' 



rA, 



+ 



tP. 



(4.6) 



(4.7) 



The equations describing the dynamics of the quasi-normal modes is obtained by elimi- 
nating A^, A^ and V^Ap from Eqs. flA.20p and flA.2ip . Equation (14. 5 p gives A^ in terms 
of A^, Ao- and P. The combination (A'^ -|- V^A^) can then be obtained, after some algebra, 
from the explicit expression of A^ and from the Hamiltonian constraint of Eq. (lA.lSp (see 
also Eq. flXSj) ). 

Thus, as discussed in appendix [Aj Eqs. (]A.23P and f lA.25P can be inserted into Eqs. 
( ]A.20p - (]A.2ip and the resulting system becomes: 



A': + 2HAi-V'A^ 

V'A, + A 



A'' + 2HA' 



A A 



A A 



.A. 



= 0, 
0. 



(4.8) 
(4.9) 
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The coefficients A^^ip, A^a and A^a = Aaip depend on the background and are: 



A, 



2d'V 1 



dip''-' 



+ =2 

Mi 



2a 



2 + 



/2 



A 
A 



A 



(Tip 



Ml 



dip \n) da 



2 + 



^' a' 



+ =2 

m; 



2a' 



,dW fo' 



da \n 



+ 2 



a 



(4.10) 
(4.11) 
(4.12) 



where, comparing with the expressions of the appendix [Aj the four-dimensional Palnck mass 
defined in Eq. (12. 5p has been introduced by trading SvrG for 1/Mp. The source terms S^p 
and Sfj appearing in Eqs. (14. 8 P and (14.90 are: 

.2 r (j2Qy. 



5. 



2MpA?^ 



2M' V-H 



cr 



?7 ^ ~^~do' 



q2 

+ y^V2(5pB-5pE), 



which are expressible in a slightly different form by using the evolution equations of y?, a 
together with the governing equation for P, i.e., respectively, Eqs. (I2.16p . (I2.17P and (13. 38 p . 
The result of this manipulation, neglecting the spatial gradients of A is: 

^2 



2MpVH 



2Mp 





[- 




Kip' 











P + Veb 
P + Veb 



^2 Q-^ 

+ y^V2((5pB-5pE), 



where Veb is defined as 



Veb = ^V2(5pB + 5pe) + V2(nB + He) - ^ ' ^'^ ^ ^\ 
3 a* 



(4.13) 
(4.14) 

(4.15) 



The system of equations derived here is the starting point for the determination of the power 
spectra of curvature perturbations to be analyzed in the forthcoming sections. 



5 Magnetized power spectra of the scalar modes 
5.1 Simplifying approximations 

The solution of Eqs. (14. 8 p and (14. 9 p during a phase of slow roll expansion determines the 
large-scale power spectra of curvature perturbations. The expressions of the coefficients of 
Eqs. (I4.10p -f l4.12p can be simplified in the limits <^ p^p and e ~ ^ 1 and it can be 
shown, for instance, that 

(X ft 
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The second inequahty appearing in Eq. (15. ip is derived by appreciating that A^a can be 
recast in the following form: 



2M? 



A. 



H 



+ 3 + ^ ^ ^ 



H^J\HJ\H 



a 



(5.2) 



and by using known identities of the slow roll dynamic^ For the simlars reasons Aaa must 
be smaller than ("H^ + Ti'). Neglecting the subleading terms in Eqs. (I4.13p - (l4.14p Sy, and 
Sa become 



where the two dimensionless variables S^n and have been introduced: 



(5.3) 



5. 



5. 



vaMp-H 



+ 



Mp dX 
X dip 
/ a' Mp dX 
\3Mpn ^ 



Mp dX\ 



5pB + 



a 



Mp dX 
SMpV. " 



5pE 



2Mpn 
2Mpn 



(5.4) 
(5.5) 



The terms containing the Ohmic current shall be neglected and the whole effect of conducting 
initial conditions will be simply encoded in the further suppression of the electric components 
as explained in Appendix [Bl 

Introducing then the rescaled variables = ax^p and = aXa and recalling that = 
V^Xip ci-iid Ao- = V^XcT) the Laplacians can be eliminated from the left and from the right- 
hand sides of Eqs. fl4.8p - fl4.9p so that the resulting equations assume the following simplified 
fori 



10. 



9. 



M 



(i'a - V q„ + =- S^{x, t) 

a Mp 



0, 



0. 



(5.6) 
(5.7) 



In the class of models introduced in Eqs. (I2.2ip - fl2.23p the slow roll parameters are given by 
e = T] = 1/a with a 3> 1 and oc a(r). Thus the coefficients A^^a and A^a become: 



A. 



2v^(3 



a 



■ifa 



A 



2(3a - 1) 



2 

1 



M 



(5.^ 



which are both suppressed in the limit e ^ 1 (i.e. a ^ 1). The coefficient A^pa is further 
suppressed because M <C Mp, as implied by the subdominant nature of a. This example 
illustrates concretely the nature of the general approximations analyzed in this section. 



^ In particular recall that (Vip/H^) = 3\/2 y^Afp and that ip/H = \/2A/p ^/e; furthermore the sub- 
dominance of a stipulates that ct <C Mp . 

""^^Note that the pump field of Eq. (|5.7p is not given by z'^/ z^r (with z^r ~ aa' /Ti). This lack of symmetry 
is ultimately related to the subdominant nature of a. 
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5.2 Primary and secondary power spectra 

In Fourier space Eqs. (15. 7p and (15. 8p become: 



z 



and their corresponding solutions are 



1 



lip 



Mp jt, 
1 



-=^S^{k,T), 



a'ir') S^{k,T')G^^\r',T)dT', 
a^iT')S^ik,T') Gi^\T',T) dr', 



(5.9) 
(5.10) 

(5.11) 
(5.12) 



where G^^\t',t) and G^^\r',T) denote the Green's function obtained from the appropri- 
ately normalized mode functions of the corresponding homogeneous equations. Denoting 
with F{k,T) and F*{k,T) the two independent solutions of the homogeneous equation, the 
corresponding Green's function is: 

F{k, t') F*{k, t) - F{k, r) F*{k, r') 



Gfc(r',r) 



W{r') 



(5.13) 



where W{t') = [F'^kr') F*{k, r') - F*'{k, T')F{k, r')] is the Wronskian of the solutions. The 
explicit form of the mode functions for and q„ are: 



Fp{k,T) 

F.{k,T) 



A/; 



^2k 

K 



2k ^ 



-kr) 



fx 



2(1 -e) 
(3-6) 



(5.14) 
(5.15) 



2(1 -e) 

The expression of the Green's function depends on the indices fi and ft of the corresponding 
Hankel functions [35l|36]. Since e ^ 1 and r] 1, the Bessel indices /i and fi can be expanded 
in powers of the slow roll parameters and ;U~3/2 + 2e + ?7 and fl = 3/2 + e. Consequently, 
to leading order in the slow roll expansion ~ /i = 3/2 and this explains why, in this limit, 
the explicit expressions of G^^\t',t) and G^^\t',t) coincide: 



— I cos Ikir' — r)l 



k'^r'r 



- + 1 sin [k{r' - r)] 



(5.16) 



Recalling that q^^\k,T) and q^\k,T) denote the solutions of the homogeneous equations 
(15. 9p and fl5.10p the primary power spectra can be computed from Eqs. fl5.14p and (I5.15P : 



{qS\k,T)q^^\p,T)) 
{q^^\k,r)q(^\p,T)) 



2n^ 
A;3 



V^{k,T)5'^^\k+p), V^{k,T) 
V.{k,T)6^^\k+ff), V.{k,T) 
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2^ 
27r2 



\F^{k,T)\\ (5.17) 
F^{k,T)\\ (5.18) 



After integration over r' the final expression can be written as: 



qaik,T) 



q^^\k,T)-Mp 



c^QB{k, t) + d^QBn{k, r) 
cMB{k, t) + d„VLBu{k, t) 



air 



(5.19) 
(5.20) 



where the sources have been evaluated to leading order in kr and 

i2B(fc, r) = =2-' ^Bn{k, r) = 



(5.21) 



The solutions fl5.19p - fl5.2ip neglects the decreasing electric modes since they are immaterial 
for the scales that had the longest time to grow and got larger than the Hubble radius 
between the last 65 and 53 efolds of inflationary expansion (see also beginning of Sec. E]). 
The rate of decrease of the electric modes is discussed in appendix[B]for the interested reader. 
The coefficients appearing in Eqs. fl5.19p and (15.201) are slowly varying functions of r 



m{f, e) 



m{f, e] 



l_ Mp/dX\ 



3Mp \V. 
3(1 -e) 



A \daj 



(l-2/)(4-2/-3e)- 



rf. = -(/,e)^Q, (5.22) 
da = m{f,e) 1 — I 



(5.23) 



A \daj' 



The function m(/, e) depends on the slow roll parameter and on the growth rate in Hubble 
units, i.e. / = J-'/'H. In the scale invariant case (i.e. / = 2, see appendix [B|) m(2, e) = 1/(3 e). 
In the pure de Sitter case and for exactly scale invariant spectrum the integration over r' 
would lead to logarithms of the conformal time coordinate which are absent in the quasi-de 
Sitter case. Equations fl5.22p and fl5.23p can be written in more explicit terms as: 



m(/,e) 
m(/, e) 



+ 1v 



3 

2e / M 



d^ = m{f, e) 7^, 



3 \Mp) ^^"K M 



da =m{f,e) (^)7<7. 



(5.24) 
(5.25) 



The results of Eqs. fl5.24p and fl5.25p assume a simple parametrization of \{ip, a), i.e. 



A((y9, a) = A* exp 



a 



^ _ 



7a + 



7^ 



'1 



2e 



(5.26) 



The second relation of Eq. f l5.26p holds, strictly speaking, in the case of power-law inflation 
where rj = e. It can be argued, however, that it remains valid in more general cases where 
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5.3 Curvature perturbations 

The curvature perturbations on comoving orthogonal hypersurfaces can be expressed in 
terms of Eqs. (l5J9|) -f l5:20|) 

as it follows from Eqs. f l4.5p -f l4T|) by recalling that = a^p' /% and z^j = aa'/Ti. The second 
equality of Eq. f l5.27p follows in the limit z^/z„ = (p'/a' ^ 1 when a is subdominant. The 
power spectrum of curvature perturbations is defined as 

(7^(fc,r)7^(p,r)) = —Vn{k,r)6^'\k + p). (5.28) 

Bearing in mind Eqs. ^26^-^28^ and (ISTTjl - dSTS]) . Eqs. (ICT]) - (15:201) can be inserted 
into Eq. fl5.27p so that the explicit expression of Vnik,T) becomes 



^2 . ^ 6B(fc,r) , QBn(fc,r) 

This expression shows that the power spectrum of curvature perturbations depends on the 
first-order correlations of the inflaton fluctuations as well as on the second-order correlations 
of the gauge fields whose related power spectra can be found in the appendix [Bl 

The first two contributions to Vnik, r) appearing in Eq. (I5.29p are the adiabatic con- 
tribution given by the inflaton and a generalized entropic contribution associated with the 
spectator field. Both C<^a(r) and V^^^r) are slowly varying functions of r (i.e. C^o-(t) — 
V^„{t) ~ 0) and are defined as 



^2 , s Mptt^fr 



' ^zI{t 



(5.30) 



What matters, for the present considerations, are those typical scales that had the longest 
time to grow and that left the Hubble radius at the onset of the inflationary phase even if, 
as we shall see, the beginning of inflation is essentially a free parameter related to the total 
number of inflationary efolds. The various contributions to the power spectrum must be 
compared the in the limit where the relevant scales are larger than the Hubble radius, i.e. 
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where, generically, )C{z) = l^^-^V^^z)!!^ so that /C(3/2) = 1 anc£| 

riad — 1 = 3 — 2/i = — 6e + 2r7, rigntr — 1 = 3 — 2/i = — 2e. (5.33) 

Barring for the dependence of the spectral index on the slow roll corrections, we can clearly 
see that Pentr(^,T) ~ e(M/Mp)^ 'Pad(fc, t). Since e < 1 and M ^ Mp the entropic contri- 
bution is strongly suppressed. If taken into account this component will lead to the kind of 
mixed initial conditions for CMB anisotropics often discussed in the literature [311 |32l |33l [31] 
also in the presence of large-scale magnetic fields [12] . 

In the simplest situation the total energy-momentum tensor of the system is conserved 
both before and after the transition between inflation and radiation [29||30]. When the stress 
tensor undergoes a finite discontinuity on a space-like hypersurface the inhomogeneities are 
matched by requiring the continuity of the induced three metric and of the extrinsic curvature 
on that hypersurface. On uniform curvature hypersurfaces the continuity of the extrinsic 
curvature is guaranteed by the continuity of A<^ and A/?. This implies also the continuity of 
v. as it can be explicitlly verified by solving the evolution equation of 7^ (see Eq. ( IA.3ip ) 
valid in the postinfiationary epoch. After the end of inflation the growth rate is zero and 
the evolution of curvature perturbations can be followed by means of a certain set of global 
variables. This discussion closely follows the considerations developed in Ref. [7]. 

In concluding this section it is appropriate to remark that Eqs. f l4.8p -( B^ and (15. 9p - 
(15.1 op . even if deduced in a specific gauge, have a gauge- invariant meaning. The gauge- 
invariant generalization of the quasi-normal modes discussed in this section is given by 

^Jf' ^ = a-Xn^ + ^^^V", (t'^ = a-Xa + Zai>. (5.34) 

Under the gauge transformation discussed prior to Eq. (14. 2 p and Xa transform as x<fi 
Xip — V^'co and X(t ^ Xa ~ '^'^o- Thanks to Eq. (14. 2 p the quantities defined in Eq. (I5.34p are 
left invariant. 

The variables of Eq. (I5.34p are the scalar field analog of the quantum excitations of an 
irrotational and relativistic fluid firstly discussed by Lukash [37] (see also [38ll39lllQ]) right 
after one of the first formulations of inflationary dynamics [41] . The canonical normal mode 
identified in Ref. [37] is invariant under infinitesimal coordinate transformations as required 
in the context of the Bardeen formalism [12] (see also [38]). The subsequent analyses of Refs. 
[13] and [H] follow the same logic of [37] but in the case of scalar field matter; the normal 
modes of Refs. [371 1131 SI] coincide with the (rescaled) curvature perturbations on comoving 
orthogonal hypersurfaces [13 US]- In the present case, as already pointed out, qj^^^ is only a 
quasi-normal mode and becomes a truly normal mode only in the case when the spectator 
component vanishes. 

-'^-'^ Recall Eqs. (|5.14p - (|5.15l) and also the well known relations among the slow roll parameters, i.e. rj = e—rj. 
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6 Growth rate of magnetic inhomogeneities 



According to the requirements spelled out in [7], the predominance of the standard adiabatic 
mode over the magnetized contributions leads to a specific bound on the magnetic field 
intensity. This logic will now be applied to the curvature perturbations induced during the 
infiationary phase with the purpose of deriving accurate constraints on the growth rate of 
magnetized inhomogeneities in Hubble units. 



6.1 Predominance of the adiabatic solution 

Demanding that the adiabatic component is dominant against the entropic and the magnetic 
contributions, Eq. fl5.29p implies 

V.Ak,r)>C^Jr)^^+V^,ir)^^, (6.1) 

having imposed M <^ Mp and, consequently, Vadik, r) ^ Pentr(^5 t)- The limit of Eq. (16.1 p 
must be applied for the scales that experienced the largest amplification. For instance the 
galactic scale crossed the Hubble radius about 53 efolds prior to the end of infiation and 
had, therefore, less time to be amplified in comparison with the scales that left the horizon 
just at the beginning of infiation. 

The bound of Eq. (16.11) is more constraining if imposed on the scales that crossed the 
Hubble radius just after the onset of infiation. Since the duration of infiation is unknown it 
is reasonable to take the total number of infiationary efolds A^t as a free parameter bounded, 
from below, by A^max denoting the maximal number of efolds that are accessible to our present 
observations (i.e. At > Ajnax)- The value of A^max is derived by fitting the event horizon of 
the infiationary phase inside the present Hubble radiuJ^ 

e-^(2..^.n,„)./^r^V'7^V"'^ (6.2) 



Ho J \H J 

where the exponent 7 controls the expansion rate during an intermediate phase ending at a 
putative scale possibly much smaller than the Hubble rate during infiation denoted by 
H. 

The parameters characterizing the dominant adiabatic component have been fixed to the 
values suggested by the best fit to the WMAP9 data |17] analyzed in terms of the vanilla 
ACDM model; this corresponds, in particular, to = 0.972 and A-ji = (2.41 ±0.10) x 10"^. 
Different data sets, like for instance the WMAP7 data |1H1 US] would imply riad = 0.963 and 
An = (2.43 ± 0.11) x 10~^. These differences are immaterial for the present considerations. 

For consistency with big-bang nucleosynthesis H,., in Eq. (16. 2 p can be, at most, 10~^^Mp 
corresponding to a reheating scale occurring just prior to the formation of the light nuclei. If 



^^For numerical estimates we recall that /iq51ro — 4.15 x 10 ^; the present value of the Hubble rate 
Ho = 100 ho Mpc"^ km/sec in Planck units is Hq ^ 1.22 x 10^^ {ho/0.7) Mp. 
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7 — 1/2 > (as it happens if 7 = 2/3 when the postinflationary background is dominated by 
dust) A^max diminishes in comparison with the case when H = Hr- Conversely if 7 — 1/2 < 
(as it happens in 7 = 1/3 when the postinflationary background is dominated by stiff 
sources) A^max increases, li Hr = H (or if 7 = 1/2) there is a sudden transition between 
the inflationary and the postinflationary regimes and, in this case, we have approximately 
A^max ^ 64 + 0.25 In e. 

In the case case of a standard postinflationary history A^max coincides with the number of 
efolds necessary to address the conventional drawbacks of the hot big bang model |50l |5T] . 
Whenever A^t > A'max the redshifted value of the inflationary event horizon exceeds the 
present value of the Hubble radius. If A^t = A^max the scales which were still larger than the 



^= 1 0"^, n =0.972 ^ = 1 0"^ n = 0.972 




Figure 1: The bound on the growth rate is illustrated by plotting the total number of efolds 
as a function of the slow roll parameter. The allowed region in the parameter space is below 
each of the various curves corresponding to different values of /. In these plots as well as in 
Figs. [2] we have set M = 10' Wp. 

Hubble radius around matter-radiation equality left the inflationary Hubble radius about 
A^max efolds prior to the end of inflation at least for a standard postinflationary history. 
Consequently the most constraining bound derivable from Eq. (16.1 p is achieved by demand- 
ing a typical number of efolds close to A^max for comoving scale of the order of denoting 
the pivot wavenumber at which the amplitude of the curvature power spectrum is commonly 
assigned when analyzing the temperature and polarization anisotropics |47l HH |49] . 

Using the results of appendix [B] and of Sec. \5\ Eq. (16. ip can be phrased as in terms of / 
(i.e. the growth rate of the magnetized inhomogeneities expressed in Hubble units). Besides 
the growth rate in Hubble units and the total number of efolds, the parameter ^ measures 
the Hubble rate in Planck units and is given by ^ = H/Mp. If the adiabatic mode is the 

^•^Some authors denote with ^ a further slow roll parameter containing four derivatives of the inflaton 
potential. The notations used here are different and, with this remark, no confusion is possible. 
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only source of inhomogeneity then there is a specific relation between e and the amplitude 
of curvature perturbations at the pivot scale qp. In the latter case and using the WMAP9 
[i7] data we 

i = ^=TiJSeAn, ^7e= (2.41 ±0.10) X 10-^ (6.3) 

IVl p 

In what follows ^ will be taken as a free parameter. The values of ^ are assigned independently 
of the values of e in all the figures of this section except for Fig. [3] holding in the case of 
Eq. (16. 3p where ^ oc y/e. The inequality (16. ip cannot be simply inverted in terms of / or in 
terms of the slow roll parameters. Equation (16. ip can instead be written as 

If^e^(^)""" >-M(/,e,r^,7.,7.,M)e^e^--(/'^), (6.4) 
where the function A^(/, e, 77, 7^, 7o-, M) is 
C-lAf, e, ^, 1^, la, M)CB{f, e)£B(/, e, Qp) +Vl^{f, e, r/, 7^, 7^, M)CBn(/, e)>CBn(/, e, gp)- (6.5) 

As already mentioned the various contributions to Eq. (16. 5p can be found in Eqs. (I5.30p . 
f lBiej) and f lB18|) - flRT9|) . 



6.2 Illustration of the constraints 

Whenever the growth rate exceeds a critical value (for a given duration of the inflationary 
A^t and for a fixed value of the other parameters), the inequalities of Eqs. (I6.ip -f l67il) are 
first saturated and then violated. With the aim of an accurate determination of the critical 
rate, the attention shall be first focussed on the case 7^ = 0; in this discussed in 

Sec. [5l the growth of the magnetic inhomogeneities is only due to the spectator field. In 
Fig. [1] the values of A"t are illustrated for different rates / as a function of e. The allowed 
region in the parameter space is below the various curves of the two plots. The vertical and 
horizontal dashed lines in the left plot of Fig. [1] correspond to a value / = 2.3 (for e ~ 0.01) 
forbidding any reasonable duration of the inflationary phase since At must be smaller than 
about 35. Larger values of / would be even more constraining for At! we conclude that the 
range of physical values is 2 < / < 2.3. For / < 2 the growth rate is not constrained by 
the predominance of the adiabatic mode since the magnetic energy density decreases (rather 
than increasing) for typical wavelengths larger than the Hubble radius. 

In Fig. [1] (plot at the right) the region of parameter space 2.05 < / < 2.15 is more 
accurately scrutinized. As the vertical and horizontal dashed lines indicate, for / ^ 2.15 
and e ^ 0.01 we are really on the borderline of the allowed region: as soon as / > 2.15, the 

^■^To avoid confusions we remind that we used throughout Alp = Afp/VSvr. Some authors prefer to use 
Mp instead of Mp and, in this case, the analog of Eq. (|6.3I) reads (H/Mp) = y/ir eAn- 
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Figure 2: The bound on the growth rate is illustrated in the planes {Nt,f) and [Nt,^). As 
in Fig. [H the allowed region in the parameter space lies below the corresponding curve. 



total number of allowed efolds drops below 60 that is insufficient to address and solve, for 
instance, the horizon problem of the conventional hot big bang model [501 EI] 

The same point is analyzed, within a complementary perspective in Fig. [2l A value 
/ = 2.2 in the left plot of Fig. |2] corresponds roughly to A^t — 30 (for e = 0.1). This means 
that the achievable number of efolds cannot be as large as Amax^ even for smaller values of 
e it turns out that At < 50. 

In the right plot of Fig. [2] the bounds on the growth rate are illustrated for / = 2.15 but 
in the plane (^, At). The line At = 60 crosses the dashed line (corresponding to e = 0.01) 
for ^ ~ 10~^Mp. Smaller values of ^ would correspond to inflationary phases occurring at 
low curvature. In this case At can be larger (for the same range of growth rates) but the 
adiabatic mode will not be able to account for the observed temperature and polarization 
anisotropies probed by direct CMB observations. The constraints set by the growth of the 
inhomogeneities are stronger than the ones simply implied by the backreaction. Given a 
growth rate / of the magnetic field we have to demand {p^ + pg) < SH"^ Mp; the latter 
condition, already discussed in Sec. [2] can be written explicitly in the case of a monotonic 
growth rate: 



doc t]c 



\h^'\x)\' + \h';^1i{x)\ 



<3H^mI, (6.6) 



327r2 

where u = 1/2 + /(1 + e) for e < 1 and where Xmin = ^minT while Xmax = ^maxT. The 
integration can be separated in the region x > 1 (where the energy density decreases as a~^) 
and the region x < 1 where the magnetic energy density increases while the electric energy 
density still decreases. By demanding that fcmin = 1/Tmin is the ffist scale leaving the Hubble 
radius at the onset of inflation and that fcmax = 1 / Tmax is the last scale leaving the Hubble 
radius at the end of inflation, the net growth of the energy density can be constrained. The 
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Figure 3: The constraints stemming from the requirement that the energy density of the 
amphfied magnetic field is subdominant in comparison with the energy density of the back- 
ground geometry. The values of ^ have been fixed as in Eq. (16. 3p . 



results are illustrated in Fig. [3] in the usual plane (A^t,e)- 

The results of Fig. |3] are less restrictive than the requirements obtained from the growth 
of the inhomogeneities. For instance, the full line of the right plot of Fig. [1] has a maximum 
for A^t — 140 while the full line of the right plot of Fig. |3] seems to allow for more than 
200 efolds. Backreaction effects can very well be under control but the inhomogeneities may 
grow larger than the adiabatic contribution. If we commit ourselves to a specific scenario, 
the magnetic energy density must not affect the equation of the spectator field. Adopting the 
parametrization previously discussed. This condition would demand 7o-/M(pb — Pe) < 3iJcr 
implying 

7a (Pb 



3ff2Mr 



Pf) M 

< ^ < 1- 

Mp 



(6.7) 



Let us consider, as an example, the case = 0. In this case 7o- is directly expressible 
in terms of / and e according to Eq. (I5.26p . The condition (16. 7p can then be plotted for 
different values of N^, e and / with the same logic leading to Figs. [T] and [H As it can 
be explicitly seen the inequalities (16.71) are satisfied with M = 1.5 x 10~^Mp for A^t — 65 
and 10~^ < e < 10~^. We therefore conclude that the predominance of the adiabatic mode 
represents a more constraining criterion than the simple backreaction requirements. 

Let us now recall that the case 7a- = 0, according to some considerations, would be 
less plausible since, in this case, the inflaton would be directly coupled to the gauge fields 
and the flatness of the potential might be in danger. In spite of these caveats, in Fig. H] 
the case 7o- = is illustrated and should be compared with Fig. [1] (obtained in the case 
7^ = 0). Provided / < 2.15 and 0.001 < e < 0.1 the total number of efolds is larger than 65. 
Conversely, larger values of the growth rate (i.e. / > 2.2) constrain the number of efolds to 
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Figure 4: The bound on the growth rate is illustrated in the plane {Nt,e). This figure is the 
analog of Fig. [T]but in the case 7cr = 0. 

be smaller than 65 and are therefore not acceptable. 
6.3 Concluding remarks 

The bounds on the growth rate can be also translated into constraints on the magnetic 
spectral index entering the phenomenological discussion of the effects of predecoupling mag- 
netic fields on the temperature and polarization anisotropics. The magnetic spectral index 
defined in [52] gives the slope of the magnetic field spectrum, i.e. Pb fc^^-i ^j^e relation 
between / and ne is ne = 5 — 2/(1 + e). The bounds on / derived in this section then 
imply a lower bound on ub > 0.6 — 4.4e. It is interesting to notice that the results of |52] are 
compatible with this limit: in a frequentistic perspective the analysis of the temperature and 
polarization correlations in the magnetized ACDM scenario implies that values < 0.9 are 
excluded to 95% confidence level. The present findings suggest that a nearly scale-invariant 
magnetic field spectrum induced by infiationary magnetogenesis is compatible with the range 
2 < / < 2.2 pinned down by requiring the predominance of the adiabatic mode during con- 
ventional infiation as argued some time ago in different contexts (see, for instance, [211 127]). 
In the nearly scale-invariant case the amplitude of the physical magnetic power spectrum 
is of the order of 1.44 x 10"^^ G (for e = 0.01 and An = 2.41 x 10"^) at the epoch of the 
gravitational collapse of the protogalaxy. 
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A Evolution equations in uniform curvature gauge 



The main set of governing equations used to derive various results discussed in the bulk of 
the paper will be illustrated in detail. With the gauge choice (14. 4p the coordinate system is 
totally fixed without the need of further conditions: because of this property the functions 
(f){x, t) and r) bear an extremely simple relation to one of the conventional sets of 
gauge-invariant variables. In the gauge (14. 4p the inhomogeneities of the energy-momentum 
tensors T^{ip), T^(cr) and Xx{p, p) are, respectively, 



5sTo^(^) = - 



4r°(a) 



Spa, 



a' 



^ + , 6,Ti{a) = -6p,6l (A.2) 



(A.3) 



where Xv aiid Xa denote, respectively, the fluctuations of the inflaton ip and of the spectator 
field a; denotes the three- velocity of the fluid in the gauge (14. 4p . The explicit expressions 
of {5p^, 6p^p) and of {6p^, 6pcr) arq^, respectively: 



/2 



,dV 



+ xV-a'|^X^),(A.4) 



SPa = 



, /2 , / / 2 



dip 
dW 



Xa (A.b) 



The fluctuations of the Einstein tensor = R'^ — R6'^/2, always in the gauge (14. 4p . are 
instead: 







1 
1 



-2{H^ + 2H') 



-H<p + {W - H^)I3 



2H(t)' 



WH<P + P' + 2HP]\51 



+ —didHp' + 2up + 



(A.6) 



(A.7) 



The combination of Eqs. f[0]) - f[0|) and f[0|l with Eqs. <^K^ - <^K1\ implies that the (00) 
and (Oz) components of the perturbed Einstein equations with mixed indices becom 



<5pt + ^Pb + 



(p + p)e + P + {p^ + p^)d^ + {p, + p,)d. 



(A.8) 
,(A.9) 



nV^fi + m^cj) = -A-nGa^ 
[W - n^)V^^ - nV^(f) = AnGa 

^^To avoid lengthy notations we wrote Sp^p and Sp^ (instead of SsP^ and 6sPa-), Sp (instead of Ssp) and 
similarly for the corresponding pressures; this notation is fully justified and unambiguous once the scalar 
nature of the fluctuations has been established, as specified by the general formulae written above. 

""^^Equations (jA.8|l and (jA.9p are commonly referred to as, respectively, the Hamiltonian and the momentum 
constraints. 



27 



where 5pt = 5p + 5p^ + 5po- and 



A, 



if' ' a 
In Eq. (lA.lOp the following practical notations 



(A.IO) 



A. 



A,, = V\., A;3 = V2/3, 



(A.ll) 



have been introduced. The {ij) component of the perturbed Einstein equations reads: 



5^ + 



+ + 2?^/3 



5pt + (^Pb + 



(A.12) 



where, in full analogy with Eq. (lA.lOp . the total pressure fluctuation 5pt has been defined: 

5pt = 6p + 6p^ + 5p„. (A. 13) 

The separation of the traceless part from the trace in Eq. (]A.12p implies the following pair 
of relations: 



{n'' + 2V.')<P + ■Hct)' + \v\(t) + 13' + 271(3) = 47rGa\6pt + 6pB + Spe) 



did^[(l) + /?' + 2np] - -V2[0 + /3' + 2'H(3]Si = 8nGa^ 



UliE) + imB) 



(A.14) 
(A.15) 

(A.16) 
(A.17) 

In Eq. ( 1A.17P the same notations established in Eq. (IA.11|) have been employed. During 
inflation the perturbative variables necessary to describe the evolution of the whole system 
are then given by A^, A/3, A^ and A^-. Neglecting the fluid sources, the Hamiltonian and 
momentum constraints of Eqs. flA.SP and (IA.9P can then be written as 



that can be further simplified by recalling Eq. f l3.25p : 

(n^ + 271') + nA'^ = A-kGo' V\6p + Sp^ + Sp^) - V' ( He + Hb 
A'p + 2HAp + A^ = 127rGa2 + Hb) . 



m'As = -ATrGa' 



V\6p^ + 5p„) + V\5pB + 6pe) 



{W - H^)Ap - HA^ = AnGaU P-^ 



ip'A^ + a' A, + {^'^ + a'')Ai, 



(A.18) 
. (A.19) 



The evolution equations of A^ and of Ao- are derived from the perturbed version of Eqs. 
(ESD and dlT]) 



A'; + 2-HA; - V^A^ + ^a'A^ + 2^a^ A, 



dip 
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i2 dX 



K + mK - V^A. + ^a^A. + 2^a^ A, 
-a'(A; + V^A,) = y |^V2(5pE - Vb). 



(A.20) 



(A.21) 



The system of Eqs. flA.20p and flA.211) can be reduced to a set of quasi-normal modes 
whose evolution equations are mutually coupled but decoupled from all other perturbation 
variables. The sum of these quasi-normal modes, weighted by coefficients that depend on 
the geometry, gives the curvature perturbations as explained in Eqs. fl4.6p -f l3?7r|) . 

The evolution equations for A^ and Ao- can be decoupled from the remaining perturbation 
variables as follows. From Eq. (lA.Op . neglecting the fluid component, we obtain an expression 
for A<^; a derivation with respect to the conformal time coordinate will give A^. The final 
result of this step is 



A,^ 


= A-nG 






Xh 


a; 


= A-nG 


\(^ 













AnGa'^ 



^VA.+ f^W+f?,)A 



n 



n 



n 



n 



2U 



p 



(A.22) 



(A.23) 



Prom the Hamiltonian constraint of Eq. (lA.Sp . always during the infiationary phase, we can 
obtain = V^A/3 and eliminate, in the derived expression, A<^ through Eq. flA.22p . This 
algebraic step leads to three typical terms: the first one contains the dependence on A^ and 
Ao-; the second term contains A'^ and A^; the third term depends on P, P' and [Sp-Q + Spe)- 
The full expression of V^A^ is: 



V^A, 
-AtiG 



AnG 
A' 



2n 

a' 
H 



A' 



+ 



AtxGo^ 
H 



dV 



A^ + 



2U + 



nJKH 



+ 



a2 dW 



n da 



2H + —jP-V'{6pB + 6pE) 



A. 



(A.24) 



By then summing up term by term Eqs. (IA.23P and (IA.24P the term (A'^ + that appears 
in Eqs. flA.20p and (IA.2ip can be explicitly obtained: 



v^As + a: 



+ 



2( 2H + 



= -AttG 
n'\ fa' 



2(2?/ + 



HJ\H 



+ 2 



.g^ dW 



A. 



P' -2^P + V\5pb + 5pe) 



(A.25) 



With the aid of Eq. (IA.25P and of the other equations derived in this appendix, Eqs. (IA.20p 
and flA.2ip reduce to Eqs. (14. 8 p and (14.90 . As anticipated these equations are mutually 
coupled but decoupled from all other perturbations variables. 
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After a transition regime, in the postinflationary phase, the couphng of the sources to 
the growth rate of the magnetic and electric fields disappears. The covariant conservation 
equation of energy-momentum tensor reduces to 



spl + (pt + pt)et + sn{5pt + 5pt) = 0, 

while the equation for the three-divergence of the total fluid velocity becomes: 



. V , [Pt+^^(Pt + Pt)].. 

^/3J H 7 : ^ [^t 

[Pt + Pt) 



A/5) 



(Pt + Pt) 



+ A„ 



0. 



(A.26) 



(A.27) 



Using Eq. (IA.26P and recalling that C = {5pt + 5pB + (5pE)/[3(pt +Pt)] the evolution of C 
becomes 



n , n 

7 : rCiPnad H : 

[Pt + Pt) Pt + Pt 



^ ) (5pB + Spe) 



P 



J-E 



3(Pt + Pt) 3(pt + Pt)a^ 



3 ' 



(A.28) 



where Spnad = ^Pt — Cgt^Pt- Similarly, the evolution equation for TZ can be almost immediately 
obtained by subtracting Eq. ( lA.SP (mutiplied by c^J from Eq. flA.14|) and by recalling the 
relation of (p to 71. The result for the evolution equation of TZ is 



7^' = + 



A7rGa'^{pt + pt)' 



(A.29) 



where Et^ is defined as 



+ 



n 



{pt + pt) {pt + pt) 



'St 



-){5pB + 5pe) + 'nE + 'nB 



(A.30) 



By taking the first derivative of Eq. (1A.29|) . the dependence on can be eliminated using 
the Hamiltonian and the momentum constraints; the final result is: by means of the other 
equations we get 



Zt 



7^ + ^ HE + nB , 



(A.31) 



where = [a^y/pt + Pt)/(^Cst). The variable ztTZ is, up to a sign, the normal mode of 
an irrotational and relativistic fluid discussed by Lukash [37] (see also EH]) with the 
difference of the source term containing the dependence on the gauge inhomogeneities. Both 
Eqs. (IA.28P and (IA.29P have been discussed in [7]. Note that, from the definition of TZ in 
terms of (f) and from the Hamiltonian constraint it turns out, as expected, that ( — TZ oc 
(see, in particular, the second paper of Ref. [12] )• So, with some caveats, the evolution of C 
can be traded from the evolution of 7Z. 
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B Second-order correlations 



The power spectra of the electric and magnetic fields measure the first-order correlation 
properties of the corresponding fiuctuations. The power spectra of the energy densities and 
of the anisotropic stresses are a measure of the second-order correlation properties of the 
electric and magnetic fields. To compute the power spectra introduced in Eqs. f l3.28p - fl3.3ip 
an explicit expression for the magnetic power spectra P-Q^q^r), P^{q,T) and Peb(q', t) is 
needed. In an exact de Sitter phase of expansion and for J-" = 21-1 = — 2/r the solution of 
Eqs. fl3.8p-f|3.1Up for the evolution of the power spectra with the correct boundary conditions 
is given bjO 

9 + 3A;V2 + A:V „ , + k^r^ ^ , fc(3 + 2A;V) 



PB(fc,r) 



47rr4 



4v]-27-2 



PEB(fc,- 



(B.l) 

The same result can be obtained directly from Eq. f l3.14p and from the related solutions in 
terms of the mode functions. In this case the wanted power spectra are: 

f d^k g3 



QB(g,r) 
QE(g,r) 
QBn(g,T) 
QEn(g, T") 



2048 vr 
2048 TT^ 
4608 7r7 
4608 7r7 



y / ^ ^[9 + 3A;V^ + fcV][9 + 3/r2 + pVlA,(g, k\ 



P 

,3 ^4 



rf3fc^(l + fcV)(l +pV)A,(g, fc), 
kp 

^-,[9 + Ser' + A;V][9 + 3pV + pV^]An(g, k), 

rC J) 



kp 



l + fcV)(l+pV)An(g,fc) 



(B.2) 
(B.3) 
(B.4) 
(B.5) 



where p = \q — k\ and the functions Ap(g, k) and Au{q, k) have been defined in Eqs. f 13.32p - 
fl3.33p . Even if the expressions of Eqs. flB.2p - flB.5p are reasonably simple, it is interesting to 
bring them to an even simpler (though approximate) form. In particular Eqs. f lB.2p -f |Bl5!) 
are equivalent to the following set of approximate expressions 



QB(g,r) 
QE(g,r) 

QBn(g,r) 
Qe 





2048 




H' 
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2048 
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4608 




H' 


i 


4608 


7r7 



Mq) ( ^ ) ^(g - -h) + - q) 



2^Bn(g) 



a* 
a 
a* 
a 



'^{q-H) + 0^{q){^ym-q) 



d{q-n) + 0^n{q)m-q) 



(B.6) 
(B.7) 
(B.8) 
(B.9) 



where the Heaviside's step function has been introduced. The factorization of the time- 
dependence has been achieved by expanding the integrands in powers of kr and pr and by 



""^^Recall that during a de Sitter stage of expansion the conformal time coordinate is negative so that ah 
the power spectra of Eq. (|B.ip are positive definite. 
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consistently keeping the leading terms in the expansion. The resulting expressions depend 
on four integrals over the momenta which can be accurately regularized and computed: 



Mq)= Jd'kkpq'A,{q,k), Oelg) = 81 J — ^A^{q,k), 

3/22 3 

IM = I d'k^j^A.iq, k), OM = I d'k^A,{q, fc), 

lBn{q)= J d'kkpq'Au{q,k), OBn(g) = 81 J^^Auiq,k) 



d^ -3 
IF 

3/22 3 

XEn(g) = / d'k^—^Auiq,k), 0^n{q) = I d^k^An{q,k). (B.IO) 
J kp J kp 

Note that Tx{q) and Ox{q) simply denote the modes of the quantity X which are, respec- 
tively, inside or outside the Hubble radius at the corresponding epoch as specified by the 
Heaviside theta functions appearing in Eqs. (]B.6|) - (]B.9p . 

The energy spectra of the electric and magnetic parts behave differently outside the 
Hubble radius. While inside the Hubble radius Qb(q', = Qe(q',t) — H^a~^, outside 
the Hubble radius Qb(q',t) — is almost constant and QE{q,T) ~ H^a~* is sharply 
decreasing. The same kind of conclusion, with slightly different numerical coefficients, also 
holds for QbhIq'^t) and QEn(9, 'T")- This means that outside the Hubble radius (which is the 
most delicate regime from the point of view of the effects on the scalar adiabatic modes), 
the magnetic components dominate against the electric ones provided the magnetic power 
spectrum is nearly scale-invariant. 

The conclusions drawn so far hold in the case of quantum mechanical initial condi- 
tions. This means that the power spectra of the electric and magnetic fields satisfy the 
corresponding equations for = — 2/r and = 0. In the case of conducting initial con- 
ditions, the situation is, in some sense, even simpler since electric fields are further sup- 
pressed at the level of the initial conditions. This means that, form the relevant equations 
of the power spectra Peb(q', t) and outside the Hubble radius O-siq, t) ~ H^a^^^^ while 
^e(q'; t) — {q / (J a^^^ where / = T jl-i] / = 2 in the case of an exactly scale invariant 
spectrum. 

The spectra of Eqs. (]B.2p - (]B.5|) are derived in the absence of slow roll corrections, i.e. 
in the case of a pure de Sitter dynamics. In the quasi-de Sitter case, the evolution equations 
of fkij) and ghij) inherit a dependence on the slow roll parameters which enter directly the 
energy spectra, slow roll corrections are then essential to derive realistic spectra and realistic 
bounds on the inflationary growth rate of the magnetic inhomogeneities. 

For typical wavelengths larger than the Hubble radius the second-order spectra including 
the slow roll corrections are given by: 

/ n \9s{<:,/) / Q \9s{<:,/) 

QB(g, r) = OB(g, /) — , QBn(A:, r) = OBn(g, /) — , (B.ll) 

f n \ 9e{^J) / a \ 9e{^J) 

QEiq, r) = O^iq, e, /) — , QEuiq, r) = OMq, /) — , (B- 12) 
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where gsie, /) and gE{e, f) are: 

gB{e, /) = 4/ - 8 + 4e /, gE^e, /) = 4/ - 12 + 2/ e. (B.13) 
The amphtudes appearing in Eqs. fIB.lip and flB.12p are: 

Oxiq, e, /) = Cxif, e) Cxif, e, q) M r"^^'^^"' (3 ^4) 

where X coincides either with the magnetic (i.e. B, BIT) or with the electric (i.e. E, EH) 
labels. In the parametrization of Eq. (IB.14P the flat spectrum of the X power spectrum 
arises for mx = 1 and the various indices corresponding to the four components are given 
by: 

msie, f) = mBn(e, /) = 9 - 4/(1 + e), mE(e, /) = mEn(e, /) = 13 - 4/(1 + e). (B.15) 
The functions Cx(/, e) are given, respectively, by: 

94/(1+^) A 

CB(/,e) = ^^^r^[/(l + e) + 1/2], CBn(/,e) = -CB(/,e), (B.16) 

24/(l+e)-2 4 

CE(/,e) = ^^^^r^[/(l + e) - 1/2], CEn(/,e) = -CE(/,e). (B.17) 

The functions Cx{f,^,q) are: 

^ ^ 8[/(l + e) + l] 8 /gy/(i+^)-4 

BlJ, e, q) 3^^^^^ + e) - 5][4 - 2/(1 + e)] 3[4 - 2/(1 + e)] [qoJ 

4 / \4/(l+^)-5 

+ ^^77^ > (B-18) 



5-4/(l + e)Vg, 



2[17-2/(l + e)] 2 2/(i+^)-4 



£Bn(/, 6, g) ^^^^^^^ + e) - 5] [4 - 2/(1 + e)] 3[4 - 2/(1 + e)] [qj 

7 f q \4/(i+e)-5 



5 - 4/(1 + e) Umax J ' ^ ^ 

^ ^ 8/(1 + 6) 8 /gy/(l+^)-6 

eU, e, gj _ 2^^^ ^ ^^^^^^^^ + e)-9] 3[6 - 2/(1 + e)] \qj 

^ ^ 2[18-/(l + 6)] 2 /gy/(l+^)-6 

EnlJ, e, gj ^^^^^^^ + e) - 9] [6 - 4/(1 + e)] 3[6 - 2/(1 + e)] Kq^J 

7 / x4/{l+e)-9 



5[9-4/(l + e)]Vg, 

The comoving scale q^ = 0.002 Mpc~^ is the usual pivot scale at which the power spectra of 
the scalar curvature are assigned. The value of go has been chosen 0.001 gp while gmax can 
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be estimated from the transition scale between inflation and radiation and it is of the order 
of 10^^ {eAnY^^ Mpc"\ The results reported in Eqs. f lRTejl - dRZT]) follow after lengthy but 



straightforward algebra from Eqs. fl3.28p - fl3.3ip . Consider, for instance, the second-order 
correlations of the magnetic energy density. From Eq. fl3.28p the explicit expression of 
2b(q', t) can be written as: 

QB(g,r) = ^^y_^rfy I„ ^s'u'\s-u\\Hi'\u)\'\Hi'\\s-u\)\'U^,s,y), (B.22) 

where y = cos is one of the angular variables arising from the integration over the comoving 
three-momentum and where the following dimesnionless vectors have been introduced 

_^ q ^ k I ^ ^1 k\ 

s = — , u=—, \s-u\ = — . (B.23) 

ah ah ah 

In Eq. flR22|) Hj;}\z) denotes the Hankel function (of generic argument z) coming from the 
solution of the mode equations including the slow roll corrections. In a specific model, such 
as the ones discussed in Sec. El will assume a specific dependence on the scale factor and 
we shall focus on the case of a monotonic dependence. The Bessel index u of Eq. flB.23P 
will then depend both on / and on the slow roll parameter. This happens since the mode 
equation for /^(r) (which is the one relevant for Eqs. (13.281) and flB.22p ) can be written aJ^ 

+ [k' - ^2 - r]h = 0, + r = a'H'lf + /(I + e)]. (B.24) 

In the present investigation we preferentially considered models, compatible with the con- 
ventional inflationary scenario, where A depends on a spectator field and it slowly increases 
during the quasi-de Sitter stage at a rate which we ought to constrain. If the slow roll 
parameters are all constant (as it happens in the case of monomial inflationary potentials, 
for instance) then aH is given by Eq. (12.200 and, to first order in e, u ~ f + 1/2 + fe. 
The integration over y in the class of integrals represented by Eq. flB.22p can be performed 
explicitly, after some algebra, when the given wavelengths are either larger or smaller than 
the Hubble radius. In connection with the lengthy algebra, Eqs. fl3.32l) - fl3.33l) imply that, 
in Eq. flB.23p . Ap{u, s, y) depends on y = cost?; the same holds for An(M, s, y) in the other 



integrals involving electric and magnetic anisotropic stresses. Using this strategy all the 
explicit expressions reported in Eqs. f lB.18p - flB.20p can be obtained after radial integration. 



^®Note that the mode function /^(t) cannot be confused with / the wavenumber has been always written 
exphcitly. With this caveat potential confusions are avoided. 
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